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| ^ ■ Abstract. We analyze the effects of lepton asymmetry on neutrino oscillations and 

on cosmological nucleosynthesis with active-sterile oscillating neutrinos. It is shown 
that small lepton asymmetries, L < 0.01, whose direct kinetic effect on nucleosyn- 
thesis is negligible, still effect nucleosynthesis considerably through their influence on 
■ oscillating neutrinos. Two different cases of lepton asymmetry are discussed: an ini- 

C*") \ tially present and a dynamically generated in oscillations. Dynamically generated in 

resonant oscillations asymmetry at small mixing angles suppresses oscillations, hence, 
the nucleosynthesis bounds on neutrino mass differences at small mixings are relaxed. 
Initially present asymmetry may suppress or enhance oscillations. The enhancement 
is a result of interchanging resonances between neutrino and antincutrino ensembles 
due to resonance waves passing through the neutrino and antineutrino spectrum. Up- 
\ dated nucleosynthesis bounds on neutrino oscillation parameters accounting for lepton 

O !~ asymmetry are presented. 

6 
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I INTRODUCTION 

The direct kinetic effect of a considerable lepton asymmetry L > 0.01 (either 
initially present or dynamically generated) on cosmological nucleosynthesis (CN) 
has been investigated in [1], and cosmological constraints on its value have been 
obtained. The effect of neutrino oscillations on CN was also studied and strin- 
gent bounds on oscillation parameters were obtained [2-7]. Neutrino oscillations, 
proceeding in the primordial plasma during CN epoch, can effect 4 He yield by 
(a) bringing additional degrees of freedom into the primordial heat bath [2], (b) 
depleting the neutrino and antineutrino number densities 1 thus slowing the weak in- 
teractions T w ~ N u El, (c) distorting the neutrino and antineutrino spectrum [3,4], 
(d) producing neutrino-antineutrino asymmetry 2 [4,8], which on its turn influences 



^ The effect was estimated for equilibrium [2] and nonequilibrium [3] oscillations and numerically 
calculated for nonequilibrium oscillations in [4]. 

2 > The asymmetry growth was estimated to be possible for great mass differences Sm 2 > I0~ 5 
eV 2 [8], however for Sm 2 < 10~ 7 eV 2 , a considerable asymmetry growth was registered in precise 
numerical studies [4]. 



the evolution of the neutrino and antineutrino ensembles and the oscillation pat- 
tern [4,9,6]. 

In this work we discuss the simultaneous effect of a lepton asymmetry and neu- 
trino oscillations on CN. Two different cases of lepton asymmetry are analyzed: an 
initially present and generated in oscillations asymmetry. We have studied the role 
of a lepton asymmetry on CN with oscillations, effective after electron neutrino 
decoupling. We have shown by a numerical analyses of the kinetics of nucleons 
and the oscillating neutrinos and antineutrinos in CN epoch, that much smaller 
asymmetries L « 0.01 exert considerable indirect effects on CN through oscilla- 
tions, due to the fact that even very small asymmetries change the medium induced 
neutrino potential energy and influence the evolution of the oscillating neutrinos. 
The prejudice that in order to influence nucleosynthesis the neutrino asymmetry 
first must grow to a considerable value L > 0.01, is not applicable for the case of 
nucleosynthesis with oscillations. 3 

Lepton asymmetry influences CN with oscillations in several ways: The neu- 
trino and antineutrino ensembles evolve differently in comparison with the case 
without asymmetry, i.e neutrino number densities, their depletion and spectrum 
distortion are changed. Also due to asymmetry term neutrino and antineutrino 
become strongly coupled and evolve differently. Lepton asymmetry changes as well 
the oscillation pattern, i.e. leads to an enhancement or suppression of oscillations. 
All these resultt into changed light elements yields, compared with the oscillatory 
case without asymmetry. This indirect asymmetry influence on CN is consider- 
able [4,5,9-11]. However, inorder to reveal it asymmetry should be considered 
selfconsistently with the neutrinos and nucleons evolution. The results obtained 
without the account of the indirect effects differ by many orders of magnitude from 
the real picture [9]. 4 

Numerical analysis, accounting for the asymmetry effect selfconsistently with 
neutrino and nucleons evolution, was provided in refs. [5,9,6,7]. See also ref. [14], 
where an precise analytical study of the asymmetry evolution, accounting for its 
back effect on oscillating neutrinos, was proposed. In the nonresonant case the 
oscillations produced asymmetry was shown to have a negligible role in CN [5]. 
However, in the resonant oscillation case the asymmetry effect on CN is consider- 
able [4,9,6,14,7]. 

In general, dynamically produced asymmetry suppresses oscillations [4,9], which 
leads to less overproduction of helium-4 in comparison with CN with oscillations 
but without an asymmetry account. Hence, the bounds on oscillation parameters 

3 ) Really in the case without oscillations the asymmetry has a direct sign-dependent effect on 
the kinetics of the nucleons, and hence, on the helium yields. Besides, when large enough it also 
contributes to the Universe energy density thus increasing the expansion rate H(t) and changing 
primordial helium-4 yield Y p [1]. 

4 ) Works considering asymmetry effect on CN in case of oscillations [12,13], providing rough 
estimates of the asymmetry growth and its effect on CN, will become more reliable when a proper 
account for the indirect asymmetry effects on CN during the full evolution of the asymmetry is 
provided in a selfconsistent analysis of the neutrino and nucleons evolution for each momentum. 



are alleviated at small mixing angles [9,6,7]. We present in sec. Ill the updated 
cosmological constraints for electron-sterile oscillation case, accounting precisely 
for the oscillations generated asymmetry. 

The effect of small initial lepton asymmetries (10~ 10 < L < 10~ 4 ) on CN with 
nonresonant active-sterile oscillations was precisely studied in [5,10,11]. It was 
found that asymmetry is able to enhance oscillations, besides its well known ability 
to suppress them [10], thus leading correspondingly to an over- or under-production 
of helium. Our analysis has shown (see sec. IV.) that the initially present asym- 
metry is able to alleviate CN bounds at large mixings and to tighten the bounds 
at small mixings (see also refs. [10,11]). 

In the next section we present the precise kinetic approach for the description of 
lepton asymmetry effects on CN in the presence of oscillations. 

II THE KINETICS 

We have used for the precise analysis the synthesis of helium-4. According to 
the standard CN the primordial helium yield depends on two compelling processes, 
determinning the nucleons freezing: Universe's cooling, H(t) ~ QeffT 2 and weak 
processes, Y w . Three neutrino flavours, zero lepton asymmetry and equilibrium 
neutrino number densities and spectrum are assumed. In the case of CN with os- 
cillations and with lepton asymmetry, all these assumptions do not work. According 
to (a)-(d): neutrino oscillations change the number of neutrino flavours; they may 
produce nonequilibrium neutrino number densities (particularly the electron neu- 
trino density may be considerably reduced in favour of the sterile neutrino density); 
the neutrino spectrum may be distorted in active-sterile oscillations; besides, even 
if initially the lepton asymmetry is assumed zero resonant oscillations may lead to 
a considerable growth of the asymmetry, which on its turn effects neutrinos and 
nucleons evolution via oscillations. This nonequilibrium picture is hard to describe 
analytically. 5 

Inorder to account for the different effects of oscillations and asymmetry on CN 
a selfconsistent numerical analysis of the kinetics of the oscillating neutrinos, the 

5 ) There exist in literature numerous analytical studies, like refs. [12,13] proposing different 
schemes for applying neutrino oscillations to solve different astrophysical or else problems in 
which schemes a central role is played by a lepton asymmetry (initially present or dynamically 
generated in oscillations) suppressing oscillations. We would like to warn that these publications 
provided too rough estimations both of the asymmetry evolution and of its effect on primordial 
nucleosynthesis. The asymmetry evolution was usually semianalitically described using differ- 
ent sorts of simplifying assumptions, in general not applicable in the nonequilibrium situation 
of active-sterile oscillations. The neutrino and antineutrino evolution was described in terms of 
particle densities (not particle density matrix) and almost equilibrium spectrum (which is not 
an adequate description of oscillationary phenomena when the growth of asymmetry is consid- 
erable). The asymmetry effect on nucleosynthesis was reduced only to L kinetic effect on CN. 
Therefore, its role in nucleosynthesis was accounted for after it has grown "enough" and hence, the 
indirect effect of asymmetry during its growth has been neglected. The asymmetry back effect on 
oscillations, in case considered at all, was assumed to be only towards suppressing oscillations. 



nucleons freeze-out and the asymmetry evolution is necessary. In our analysis the 
set of the following coupled integro-differential equations describing the evolution 
of the neutrino p and neutron number densities n n was solved simultaneously and 
selfconsistently: 

Mt) _ R dpit) | 
dt v dp v 

+ % [H , pit)} + iV2G F (C - iV 7 [ a , p(t)} + O (G 2 F ) , 

dm _ d-pjt) | 

dt v dp v 

+ % [n o , pit)} + iV2G F (-£ - Q/Ml) iV 7 [a, p(t)} + O (G 2 F ) , (1) 
(dn n /dt) = Hp n (dn n /dp n ) + 

+ J dil(e~,p, is)\A(e~p -> un)\ 2 [n e -n p (l - p LL ) - n n p LL (l - n e -)\ 

-J dfi(e + ,p, v)\A{e + n ^ pv)\ 2 [n e +n n (l - p LL ) - n p p LL (l - n e +)\ . (2) 

where = U* e Uj e , mixing just in the electron sector was assumed z/j = UuVi{l = 
e, s). p v is the momentum of electron neutrino, n stands for the number density of 
the interacting particles, dQ(i,j, k) is a phase space factor and A is the amplitude 
of the corresponding process. 

These equations provide simultaneous account of the different competing pro- 
cesses, namely: neutrino oscillations, Hubble expansion and weak interaction pro- 
cesses. 7i is the free neutrino Hamiltonian. The 'nonlocal' term Q arises as an 
W/Z propagator effect, Q ~ E v T '. £ is proportional to the fermion asymme- 
try of the plasma and is essentially expressed through the neutrino asymmetries 
£ ~ 2L Ve +L Vli +L„ T , where L^ T ~ {N^-N^/N^ and L Ve ~ J d 3 p(p LL -p LL ) /iV 7 . 

The neutron and proton number densities, used in the kinetic equations for neu- 
trinos, are substituted from the numerical calculations of eq. (2). On the other 
hand, pn and Pll at each integration step of eq. (2) are taken from the simulta- 
neously performed integration of the set of equations (1). 

The equations are for the neutrino and neutron number densities in momentum 
space. This allows to account precisely for the spectrum distortion effect and 
neutrino depletion effects of oscillations, as well as to follow the evolution of the 
neutrino asymmetry and its back effect at each neutrino momentum. 

In our numerical analysis the spectrum distortion was described by 1000 bins for 
the nonresonant case and by 5000 bins for the resonant case. In case the spectrum 
was described by N bins, a system of 6iV + 1 coupled integro-differential equations 
following from (1) and (2), was numerically solved. 

The numerical analysis was provided for the characteristic temperature interval 
[2 MeV, 0.3 MeV] and the full set of oscillation parameters of the active-sterile 



Figure 1. Cosmological constraints 
for the electron-sterile neutrino os- 
cillations, are presented by the solid 
curves Y p = 0.24. The dashed curve 
shows the contour without asym- 
metry account. The dotted curve 
shows solar neutrino LOW solution. 



Figure 2. The dependence of 
the neutron number density rela- 
tive to nuclcons X n = N„/(N P + 
N„) for the case of oscillations with 
sin 2 (2i?) = ltr - 05 on the value of 
the initial asymmetry is plotted [10]. 



Figure 3. On the 5m? — -d 

plane the isohelium contours Y v = 
0.24 calculated in the discussed 
model of CN with neutrino oscil- 
lations and initial lepton asymme- 
tries L = 10~ 6 (solid curve) and 

L = 10~ are shown. The dot- 
ted curve denotes LOW solution. 



oscillation model [4]. We calculated precisely the n/p-freezing, essential for the 
production of helium, till temperature 0.3 MeV, and accounted adiabatically for 
the following decays of neutrons till the start of nuclear reactions at about 0.1 MeV. 



Ill DYNAMICALLY GENERATED LEPTON 

ASYMMETRY 

Our numerical analysis showed that in the resonant oscillation case the dynamical 
neutrino-antineutrino asymmetry grows up to 4 orders of magnitude. I.e. starting 
with asymmetries of the order of the baryon one it reaches maximum a value 10~ 5 . 
Hence, having in mind this small value the registered asymmetry effect is totally due 
to its indirect influence on CN via oscillations. Dynamically produced asymmetry 
at small mixing angles suppresses oscillations, which leads to less overproduction of 
4 He in comparison with CN with oscillations but without an asymmetry account. 
Hence, the cosmological constraints on oscillation parameters are alleviated at small 
mixing angles. 

The updated constraints on active-sterile neutrino oscillations, precisely account- 
ing for the asymmetry generation, spectrum distortion and the depletion of the 
neutrinos, are presented in fig. 1. The plots correspond to Y p = 0.24. The net 
indirect asymmetry effect on CN is given for the resonant case by the difference 
between the dashed curve (without asymmetry account) and the solid one. Due 
to asymmetry growth account, and the corresponding suppression of oscillations, 
Y p overproduction is not so strongly expressed at small mixing angles, hence CN 
constraints are alleviated at small mixing angles for 5m 2 < 0. 

In the resonant case the cosmological constraints at large mixings are \5m 2 \ < 
8.2 x 10~ 10 eV 2 . In the nonresonant case 5m 2 > an analytical fit to the exact 
constraints is: <5m 2 (sin 2 2 r d) A < 1.5 x 10~ 9 eV 2 . The constraints in both cases are 



strengthened compared to the previous ones [2] due to the precise account of the 
spectrum distortion and to the exact kinetic approach to the neutrinos and nucleons 
evolution. 

According to these constraints the LOW active-sterile solution to the solar 
neutrino problem, which is favoured by the analysis of the recent experimen- 
tal data of total measured rates and day and night spectrum measured by Su- 
per Kamiokande [15] is almost completely excluded [11]. 

IV INITIALLY PRESENT LEPTON ASYMMETRY 

The role of initially present relic asymmetries on CN with nonresonant active- 
sterile oscillations was precisely studied, following the lines of work described in 
sec. II. A wide range of L values [1CT 10 — 1CT 4 ] was analyzed. Such small asym- 
metries have only indirect effects on CN. On fig.2 the dependence of the produced 
helium, in a CN model with oscillations on the initial asymmetry is plotted. It 
was found that the asymmetry is able also to enhance oscillations, besides its well 
known ability to suppress them. The enhancement is a synthetic effect of a resonant 
wave passing through neutrino spectrum till lepton asymmetry changes sign and 
followed by a similar 'spectrum' resonance at the antineutrino ensemble and vice 
versa [10]. Hence, this enhancement has a complex spectral character and could 
be revealed only by a precise kinetic approach to the oscillation problem. The 
influence of the asymmetry should not be reduced only to oscillations suppression, 
as usually believed. Depending on the concrete values of oscillations parameters 
the asymmetry may suppress, enhance or not influence oscillations, thus leading 
correspondingly to an under- , over-production of helium or not change its abun- 
dance. Hence, in order to judge the real asymmetry effect on CN, precise numerical 
analysis for the concrete oscillation parameters and L values is obligatory. Initially 
present asymmetry may relax CN bounds at large mixings and tighten the bounds 
at small mixings. 

Qualitatively, for oscillations affective after the freeze-out of the electron neu- 
trino, the asymmetry effect is as follows: L < 1CT 7 has negligible effect on CN; 
1CT 7 < L < 1CT 5 enhances oscillations due to the spectrum wave resonance [10], 
resulting into an enhanced overproduction of helium-4; L > 10~ 5 leads to a sup- 
pression of oscillations and relaxation of the CN bounds (fig. 2). 

In fig. 3 the isohelium contours Y p = 0.24 are presented, for different L: L = 
10~ 6 and L = 10~ 10 . For small mixing angles the asymmetry L = 10~ 6 enhances 
oscillations, which reflects into stronger bounds on oscillation parameters, while for 
large mixings this asymmetry suppresses oscillations and CN bounds are weakened 
compared with L = 10~ 10 case. 

A similar investigation, of initial asymmetry effect on CN with oscillations, for 
the resonant case will be a more complicated task due to technical problems: The 
usual explicit numerical approach is not applicable for the description of the asym- 
metry evolution, because the neutrino evolution equations at resonance have high 



stiffness; besides the resonance case deserve much greater number of bins for the 
spectrum distortion description. To solve the stiff equations numerically, implicit 
methods should be used. For 5000 bins of the spectrum a system of 30000 equations 
describing the neutrino evolution should be solved simultaneously. However, this 
investigation is interesting, as far as such small values of the initial relic asymme- 
try are not excluded neither from observations nor from some profound theoretical 
principle. 

In conclusion we would like to stress that small asymmetries, initially present or 
dynamically generated, influence CN thanks to their backfeed effect on oscillating 
neutrinos. The account of spectrum distortion of the oscillating neutrinos as well 
as the selfconsistent account of neutrinos and nucleons evolution is essential for re- 
vealing the indirect effect of small lepton asymmetries. Therefore, a precise kinetic 
approach should be provided when analyzing the effect of lepton asymmetry on CN 
with oscillations. 
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atmosphere of the conference and for the financial support of the participation of 
D.K. M.C. thanks Theory Division of CERN where this work was prepared. 
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